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Abstract The most usual formulation of the Laws of Thermodynamics turns out to be suitable for 
local or simple materials, while for non-local systems there are two different ways: either modify this 
usual formulation by introducing suitable extra fluxes or express the Laws of Thermodynamics in terms 
of internal powers directly, as we propose in this paper. 

The first choice is subject to the criticism that the vector fluxes must be introduced a posteriori in 
order to obtain the compatibility with the Laws of Thermodynamics. On the contrary, the formulation 
in terms of internal powers is more general , because it is a priori defined on the basis of the constitutive 
equations. Besides it allows to highlight, without ambiguity, the contribution of the internal powers in 
the variation of the thermodynamic potentials. 

Finally, in this paper, we consider some examples of non-local materials and derive the proper 
expressions of their internal powers from the power balance laws. 
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1 Introduction 

The modern approach to the thermodynamics of continuous systems, known as Rational Thermody- 
namics, traces back to the works of Coleman, Gurtin, Noll and Truesdell [B], [IS], [30], [33] . It's our 
opinion that this important new approach in the study of thermodynamic phenomena has achieved its 
formal stage in the papers of Noll [3Tj and Coleman and Owen [7J. As it is well known, this innova- 
tive point of view distinguishes the field equations, which hold for all systems, from the constitutive 
equations, which characterize the materials. Moreover, according to this school, restrictions on the 
constitutive equations can be obtained from the Principles of Thermodynamics. 

In most of the works of the previously cited authors only simple materials are considered. These 
systems are defined by Truesdell ([33], pag.35) as those materials for which everything there is to know 
can be found out by performing experiments on homogeneous motions of a body of those materials, from 
whatever state they happen to find in it. In the same book, Truesdell gives also a formal definition of 
simple materials, in presence of purely mechanical and isothermal processes, through the constitutive 
equation T(x,t) = T(F*(a;, ■);%) which links the stress tensor T to the history up to time t of the 
strain gradient F, at any point x of the body B. 

A more rigorous and general definition of simple material can be found in |31) and [JJ. These 
materials are always characterized by the fact that the First Law of Thermodynamics is given by 

p^e = T.L- V-q + pr, (1.1) 
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where p is the mass density, e the internal energy, L the velocity gradient, q the heat flux and r the 
heat supply, while the Second Law of Thermodynamics is formulated through the Clausius-Duhem 
inequality 

P3T?7>-V- ^ (1.2) 



dt 

•q being the specific entropy density and 6 the absolute temperature. 

The expressions (jl.lj) and (|1.2j) of the First and Second Law of Thermodynamics cannot be applied 
to those systems, called non-local or non-simple materials, for which the constitutive equations for T 
and q depend on the gradients of order higher than the first of the thermo- kinetic variables (F, 9). 

In order to understand how it is possible to give a formulation of the First Law of Thermodynamics 
for a wider class of materials, we consider the heat balance equation [7] 

ph = -W-q + pr, (1.3) 

where h is the rate at which the heat is absorbed per unit mass. By comparing (jl.3p with (|1.1[) . the 
First Law of Thermodynamics for simple materials can be rewritten as follows 

P j t e = T-L + ph. (1.4) 

Moreover, following many classical texts of Thermodynamics (Truesdell [35], Gurtin [17] . Germain [TBI 
116] , etc.), we introduce the internal mechanical power and the internal heat power V l h = ph. Then 
the balance equation for the energy (|1.4[) assumes the classical form 

Pj t e = V: n +Vl (1.5) 

which is more general and can also be applied to non-simple materials. Obviously it coincides with 
(li~4|) when 

K = T ■ L, (1.6) 

but this choice, which is surely correct for simple materials, cannot be assumed in general. For this rea- 
son Dunn and Serrin [5] suggested to change the formulation (jl.ip of the First Law of Thermodynamics 
with the new expressiorQ 

p^-e = T- L- V- N- V- q + pr, (1.7) 

dt 

where N is called interstitial work flux. Expression (11.71) has been actually introduced by the authors for 
some non-simple materials including the elastic materials of grade three and the subclass of materials 
of Korteweg type, but its framework is more general ( [13li27| ) . 

In this paper we start from the statement that the representations (|l.ip or (ll.4[) of the First Law of 
Thermodynamics cannot have universal character and propose the more general classical formulation 
(|1.5[) . where the internal mechanical power will assume different expressions depending on the 
materials under study. 

It is also interesting to observe that a comparison between Dunn and Serrin's formulation (|1.7[) and 
the classical one (jl.5p gives 

^ = T-L-V-N, (1.8) 

which turns out to be an indirect way to define V % m . Let us consider, as an example, the constitutive 
equation for the stress tensor of a second grade material, i.e. 



= J-2 



V-T, 



where T2 and T3 are respectively second and third order tensors. The internal mechanical power is 
given by ( ) 

V l m = T 2 • Vv - T 3 • V(Vv). (1.9) 

Since 

T 3 • V(Vv) = V • (T 3 Vv) - (V • T 3 ) • Vv, 

1 Only recently it has been observed that expression (|1.7I) and the introduction o f th e extra-flux N in the 
First Law of Thermodynamics had been already proposed in 1982 by T. Manacorda [25], see also |26] . 
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expression (|1.8p turns out to be equivalent to (|1.9p introducing the vector 

N = T 3 Vv. 

Nevertheless, we prefer the representation of the internal mechanical power (jl.9p . since it is in agreement 
with the expression of the internal power as a differential form and since (|1.8[) results a " hybrid" form 
because it contains a term of ambiguous differential type, T • L, and a flux term. 

Similar arguments can be applied to many non-local mechanical materials, to the study of the 
phase transitions of superconductors and superfluids, to the transition of a binary mixture satisfying 
the Cahn-Hilliard equations and, in electromagnetism, to dielectric media with quadrupoles. 

In literature some authors maintain the expression for the First Law of Thermodynamics 

and, in presence of the above non-simple materials, modify the Clausius-Duhem inequality (jl.2p . They 
consider, as suggested in 1967 by Miiller [35], the following inequality 

Pj t r]>-V-3> + p r -, (1.10) 
where the vector <fr, still called entropy flux, can be defined by 

* = *' + §, (1-11) 

being the entropy extra-flux. 

This approach is useful only in presence of non-local effects on the heat flux. Instead it is not fruitful 
in the study of materials with constitutive equations describing non-local behaviors in mechanics, 
electromagnetism, chemistry, etc.. For these materials the Second Law assumes the classical form 
(|1.2p . while we need to provide a suitable expression for the (mechanical, chemical, electromagnetic, 
etc.) internal powers to use in the First Law. 

On the other hand, some thermal effects are well described by non-local models for which the 
Clausius-Duhem inequality (jl.2D must be substituted by (jl.lOp . An example is the Guyer-Krumhansl 
equation [2TJI|2"T] . 

4q + - = -k(0)V0 + a(V ■ Vq + 2VV ■ q), (1.12) 

at t 

which has beeen derived by solving the linearized Boltzmann equation of phonon gas hydrodynamics. 
Here r, k and a are suitable constitutive coefficients and this model, as proved in [3], satisfies the 
inequality (jl.101) by choosing 

*' = -((Vq)q + 2qV-q). 

K 

The presence of the extra- fluxes N in (|1.8p and €>' in (|1.10p - (jl.lip opens to the criticism of a 
formulation of the Principles of Thermodynamics which turns out to be not universal, since these 
vectors must be defined a posteriori. 

As we have already observed, a general formulation for the First Law of Thermodynamics is given by 
(|1.5p , where the internal power densities V % m and V % h are a priori defined on the basis of the constitutive 
equations. In order to get a similar general expression for the Second Law of Thermodynamics, we 
introduce the entropy balance equation 

^ = ^ V - q + 4 (L13) 

which can be derived multiplying equation (|1.3[) by the coldness 6 . 

This equation distinguishes between the contributions of the internal and external actions because 
(|1.13p can be rewritten as 

( * + £ q . v »__v.(s) + ,j. (i.i4) 

For simple materials in the heat flux (i.e. for materials with a local constitutive equation for the heat 
flux), we define internal entropy action the left-hand side of (j!.14[) . which is a differential form, and 
external entropy action the right-hand side of (|1.14p . i.e. 



(1.15) 
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Therefore, thanks to (jl.141) . the usual Clausius-Duhem inequality (|1.2p assumes the equivalent form 

(1.16) 



d h 1 „„ 



or, by virtue of (|1.15[) . it can be expressed only by means of the internal entropy action as follows 



(1.17) 



We observe that the expression (j 1 . 16|) holds only for simple materials in the heat flux, while the relation 
(|1.17|) is consistent also with the modified Clausius Duhem inequality (ll.lOp when the internal entropy 
action is given by 

h 1 



A' 



r q- ve- V- 



(1.18) 



It is therefore clear that it is necessary to give a proper definition of internal and external entropy 

action. In fact, if for example we consider a heat conductor of Guyer-Krumhansl type (I1.12[) with 

c 

k(9) = — , the entropy balance equation becomes 
9 



-V 



q 


l 




+ -q- 


.9. 


c 



4-q + - - a(V • Vq + 2VV • q) 

at t 



Since 

we obtain 

h d 



q • (V • Vq + 2VV • q) = V • [(Vq)q + 2qV • q] - | Vq| 2 - 2| V ■ q| : 



P dt V 2c 



|:)-^ 2 -^(|v q | 2 + 2|v.q| : 



§ + -((Vq)q + 2qV-q) 

V c 



r 

Pa- 



It follows that, for such material, the internal entropy action is defined by 



A' 



Pi 



1 



dt \ 2c j ct 
while the external entropy action is given by 



|Vq| 2 + 2|V • q| : 



((Vq)q + 2qV-q) 



(1.19) 



Finally, in this paper, we suggest to describe the evolution of a thermodynamic system through the 
virtual mechanical power and entropy action. To this end, it is important to introduce the domain 
Q C M 3 on which are defined the fields and the processes of the virtual velocity v and of the virtual 
temperature 9. By denoting the internal and external virtual mechanical powers respectively with 
Vln, and the internal and external virtual entropy actions with A l en , A e enl we are now able to write 
the balance equation for the virtual mechanical power 



p— v • v dx 

dt 



VLdx= \ V^dx 



and the balance equation for the virtual entropy action 



In presence of simple materials, the previous equations become 



/ pjv-vdi+ / T-Vvdx= [ [V • (Tv) + pf • v\ dx, 



ve 


dx — I 


-V • 


(f) 


r 

+ % 


dx 




Jn 











where f denotes the external forces and the symmetry the the stress tensor T is assumed. 
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2 Definition of thermo-mechanical material as dynamic system 

The thermo-mechanical properties of a material are based on the notion of state and process. Here we 
provide these notions following Noll and Coleman and Owen [7] ( see also 

We consider a body occupying the placement £>. For any material point l£6,a thermo-mechanical 
process P in ^"of duration dp > 0, is a piecewise continuous function on [0, dp) with values in E, open 
and connected subset of a vectorial space. 

The notation Pft^ta) denotes the restriction of P to [ti, ta) C [0, d p ). In particular, we denote by P t 
the restriction of P to the interval [0, i), t < dp. Let Pi and P2 be two processes of duration dp 1 and 
dp 2 , the composition Pi * P2 of P2 with P± is defined as 



p + p - J p iW' *e[0,diO 
12 \P 2 (t-dp 1 ),te [d^, dp, +dp 2 ) 



We can therefore give the following definition of thermo-mechanical material as dynamic system. 

Definition 2.1 Let Bt the placement occupied by the body at time t. A thermo-mechanical system, 
at any t and x £ B t , is a set {77, Vt, S, % U} such that 

1. 77 is the space of the thermo-mechanical processes P satisfying the following properties: 

i ) if P £ 77, then P[ tl ,t a ) G 77 for every [ti, t 2 ) C [0, dp), 

ii ) if Pi, P 2 £ 77, then Pi * P 2 £ 77, 

2. Vt is the space of the virtual thermo-mechanical processes Vt, i.e. the space of the processes that 
may be specified independently of the actual evolution of the system at some arbitrarily chosen but 
fixed time t, 

3. S is an abstract set, called state space, whose elements a are called states, 

4. f :£x/I- > Z 1 is a map, called state transition function which assigns, to each state a and process 
P, the state at the end of the process and satisfies the following relation 

XO, Pi * P 2 ) = x{x{<J, Pi), P 2 ) , a £ S , Pi, P 2 £ 77, 

5. for any state cr, process P of duration dp, time t £ [0, dp) and virtual process Vt £ Vt, it is defined 
the response function at time t 

U t : £ x E xV t 4lxl 

which assigns the internal virtual mechanical power V % m and the internal virtual entropy action„4* n , 
namely 

U(t) = Ut(a t ,P(t),V t ) = (V: n ,Al n ) t t £ [0,d P ), 

where a t = x( cr , Pt)- 

Remark 2.1 Because of the linear dependence of the virtual power and action on the virtual processes, 
the function U(t) is linear with respect to v t . 

Remark 2.2 The knowledge of V l m and A\ n allows to define the internal mechanical power V % m and 
the internal entropy action A\ n . Moreover, if we restrict A l en to spatial homogeneous thermal processes, 

we can also define h and the heat flux power V q — A l en — p—. 



3 Thermodynamic Laws 

In order to introduce the Thermodynamic Laws, following [7], we introduce the notion of cycle. 
Definition 3.1 A pair (a, P) £ S x 77 is called a cycle if x(<7, P) = a. 

First Law of Thermodynamics. In any cycle (a, P) the sum of the rate at which the heat is absorbed 
by the body and the power of the internal forces vanishes, i.e. 

j> " {ph(at,P(t))+P l m (a t ,P(t))} dt = 0. 
A consequence of the First Law of Thermodynamics is the following theorem 
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Theorem 3.1 The First Law of Thermodynamics implies the existence of a unique (up to an additive 
constant) continuous function e : S — > M, called internal energy, such that, for any pair of states ax, 
(72 and process P with x(°~l > P) — °~2 , we have 

P e(a 2 ) - pe(ax) = f " {ph(a t ,P(t)) + V l m (a t , P(t))} dt. 
Jo 

Moreover, under suitable hypotheses of regularity on the constitutive functionals, 

P j t e{a t ) = P h(a u P(t)) + V l m {o t ,P{t)) (3.1) 
for any instant t of continuity of the right-hand term in (|3.1[) . 

Second Law of Thermodynamics. In any cycle (a, P) the internal entropy action is non-positive, 
i.e. 

(■dp 

Al n (o- t ,P(t))dt < 0. 

As a consequence of the Second Law of Thermodynamics we get the following theorem 

Theorem 3.2 The Second Law of Thermodynamics implies the existence of at least one upper semi- 
continuous function 77 : S — > K, called entropy, such that, for any pair of states o~\ , 02 and process P 
with x(<7i, P) — a~2, we have 

pr}(o- 2 )-pv(o-i)> [ Al n (a t ,P(t))dt. 



Moreover, under suitable hypotheses of regularity on the constitutive functionals, 

pf t v(a t )>AUa t ,P(t)) (3.2) 



for any instant t of continuity of the right-hand term in 

When only isothermal processes are involved, the two Laws of Thermodynamics are replaced by the 
Dissipation Principle of mechanical work. 

Dissipation Principle. In any cycle (<r, P) the internal mechanical power is non-negative, i.e. 

dp 

V l m {a u P(t))dt>Q. 
A consequence of the Dissipation Principle is the following theorem 

Theorem 3.3 The Dissipation Principle implies the existence of at least one lower semicontinuous 
function ip : S — > K, called free energy, such that, for any pair of states o~\, &i and process P with 
x(ci,P) = 02, we have 

p^{u 2 )-(4{o- x ) < [ " V l m {a t ,P{t))dt. 



Moreover, under suitable hypotheses of regularity on the constitutive functionals, 

Pj^{a t )<P l m {a u P{t)) (3.3) 
for any instant t of continuity of the right-hand term in 
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4 Virtual powers and thermo-mechanical materials 

For a thermo-mechanical system, we define the field of the generalized virtual motions at time t as the 
pair 

1 \ 



?*(•)= v t (-), 



1 



where vj(-) and — — denote respectively the vector field of the virtual velocities and the virtual 
#*(•) 

coldness on B t . 

The virtual powers of the internal forces at any x in Bt are taken as linear forms in the gradients 
of the virtual motions v t evaluated at x. For a large class of materials we can represent the internal 
virtual mechanical power as 

P* m (t, x) = J2 T l+ i (t, x) ■ V«v t (z) (4.1) 

i>l 

where V*-™- 1 denotes the gradient of order n and Tfc is a tensor of order k. Moreover, the internal virtual 
entropy action is defined by 



A\ n {t, x) = rJ—g(t, x)h(t, x)-Y\ x ) ■ V 
Ot{x) 



(i) 



e t {x) 



(4.2) 



where is a tensor of order i. Finally, the virtual mechanical power and entropy action of the external 
and contact forces are respectively given by 



r^(t,x) = v 



■A en [t, x) 



i>i 
\i>l 



+ p(t,x)f(t,x) ■ v t (x), 



(x) 



+ p{t, x)r(t, x)-. 



t(x) 



The principle of virtual power is the requirement that the external and internal virtual powers and 
entropy actions be balanced, i.e. 



Q (t,x)-v(t,x)-^ t (x)+V l m (t,x)=V^(t,x), Ai n {t,x)=A e en (t,x) (4.3) 

for any generalized virtual motion. 

From now on, whenever no ambiguity arises, the dependence on x and t will be understood and 
not written. 

Definition 4.1 A thermo-mechanical material is called a simple material if 

T 4+1 = , q, ; = , i > 2 
for any time t and x G Bt ■ Otherwise it is called non-simple material. 

Remark 4.1 From the previous definition, it follows that, in presence of a simple material, (|4.ip ~ (|4.2p 

become 



£t(<r,P,7t)=T 2 (a,P)-Vv t , 
A\ n (a, P, 7t) = ph(a, P) I - q : (a, P) ■ V 



where T 2 is the classical Cauchy stress tensor and qi is the Fourier heat flux vector. 
It should be noticed that the arbitrariness of the virtual process yields that the triplet (T2,h, qi) 
is uniquely determined when the internal virtual power and entropy action are known. This is in 
agreement with the well known rational thermodynamic theories ([S], [B], [7], [30], [SI], [33], [33], [33]). 
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By using the classical notations (T, h, q) instead of (T2,/i,qi) and (L,#, g) instead of (Vv, 9,V9) 
(where the dot denotes the time derivative), the relations (|3.1j) and p.2j) become 

^efo) = h(a t ,P(t)) + -T(a t ,P(t)) ■ L(t), (4.4) 
ere p 

>^TT + ?(i)*' P( ^* (4 - 5) 

Moreover, introducing the free energy potential 

if) = e-9r), (4.6) 

P~4]) and (@3]) yield 

|^K) < -T(a t ,P(t)) ■ L(t) - -Lq^t, P(t)) • g(t) - 0(f)»j(0t). 

In the context of the theory of simple materials it is possible to obtain mathematical models well 
describing the thermo-mechanical behavior of several materials. Here we recall the most elementary 
examples. 

— Classical Thermoelastic Systems 

The state of these materials is given by the deformation gradient F and the absolute temperature 
9, i.e. a — (F,9), and the constitutive equations are: 

T = T(ct), h = A(a) - L + B{a)9, q = K(a)g. 

— Newtonian Viscous Fluids 

The state of these materials is given by the mass density g and the absolute temperature 9, i.e. 
a = (g,6), and the constitutive equations are: 

T = [—p{a) + A(er)trD]I + 2p(cr)D , h = Ai{a)g + B x {a)8 , q = K (cr)g, 

where D = sym{L}. 

— Materials with Fading Memory 

The state of these materials is given by the histories of the deformation gradient F* and of the 
absolute temperature 9 t , i.e. a — (F',0*) , and the constitutive equations are: 

T = T'{a), h = A'(cr) • L + B'{a)9 , q = K'(a)g. 

— Materials with internal variables 

The state of these materials is given by the deformation gradient F, the absolute temperature 9 
and the internal variables a 6 K™, i.e. a = (F,9,a), and the constitutive equations are: 

T = T*(a), h = A* (a) ■ L + B*(a)9, q = K*(a)g. 



5 Non-Simple Materials 

As already observed in the Introduction, a significant difference between simple and non-simple mate- 
rials is pointed out when the mechanic power and entropy action are defined. For simple materials, the 
internal mechanical power V % m is always expressed by (jl.6l) and, therefore, the First Law assumes the 
form or (|1.4[) . while the internal entropy action A l en has always the form (|1.15[) and the Second 
Law is given by the Clausius Duhem inequality (|1.2p . 

For more general materials there are two ways: either (jl.ip and (|1.2[) must be modified introducing 
suitable extra fluxes, as proposed by Dunn and Serrin [8] and Miiller [29 , or the expressions of the 
internal mechanical powers and of the internal entropy actions V % m and J¥ en must be properly derived 
from the power balance laws. 

If we choose the first option, the First and Second Law of Thermodynamics are respectively ex- 
pressed by (|1.7p and (jl.lOp . It should be noted that the internal powers and actions defined in (|1.8p 
and (|1.18p are not given by a differential form, but by a hybrid representation containing also external 
fluxes. In any case, we prefer the second option also because it allows us to highlight, without am- 
biguity, the contribution of the internal mechanical power and entropy action in the variation of the 
thermodynamic potentials. 
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5.1 Second grade materials 

In this subsection we limit our attention to second grade materials. 
Definition 5.1 A thermo-mechanical material is called a second grade material if 

T i+1 = , q, = , i > 3 

for any time t and x G Bt- 

In the sequel we will show how through our theory it is possible to obtain a framework in agreement 
with the thermodynamic principles. 

5.1.1 First Law of Thermodynamics 

For second grade materials, taking in account Definition 15. 1[ (|4.1[) becomes 

K>, P,ut) = T 2 (a, P) ■ Vv t + T 3 (a, P) ■ V^v t , 

This expression fully agrees with the wide literature on these materials. Among all the papers, we 
recall the one by Fried and Gurtin [13]. Here the authors, in accordance with [16], introduce a general 
virtual power principle by asserting that the internal power has the form 

W mt = J (T 2 ■ Vv + T 3 • V^v) dx, (5.1) 

where T 2 and T3 are respectively second and third-order tensor£0. Moreover they observe that it is 
reasonable to modify also the classical kinetic energy and introduce a generalized kinetic energy 

K* = i( P v 2 +/3|Vv| 2 ) 

to take into account its dependence on the gradient of velocity. 

This last observation allows us to study, within this theory, also the Kirchhoff thermoelastic plate. 
Following [TU], [22], [23], the constitutive equation for the stress tensor is 

T = -a\7 (Au) + cV9 + T K (5.2) 

where A denotes the Laplace operator, u is the vertical displacement, Tk = bVil is called kinetic 
stress and a, b and c are suitable material constants. 
Multiplying the equation of motion 

pu = V • T + pf 

by u, we obtain 

Jt Q^ 2 ) =V-(Tu)-T.Vti + p/« (5.3) 

and making use of the constitutive relation f|5.2[) the balance equation of the mechanical powers assumes 
the following form 



ld_ 

Ydt 



pii 2 + a (An) + b (Vu)' 



- c£ A u 

= V • {[-oV (Am) + bWil + cV6] u+(aAu~ 0)Vu} + pfu. (5.4) 
By virtue of (|4.3|) i the internal and external mechanical powers are therefore given by 



1 d 

V = 

m 2 dt I 



a(Auy + b(Vii) -cOAii, V m = -V • N' + pfu , (5.5) 



2 The representation (|5.1[) can be applied to a solid or liquid flow at small length scale and, with minor 
modifications, also to second grade nanocristaline elastic materials [14| . 
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where 

N' = [oV (Am) - bVil - cV0] u - (a A u - 6) Vii. 
Consequently, the local form (|3.ip of the First Law becomes 



d Id 

P dt 6 ~ 2dt 



a(Auf + b(Vuf 



c8 Aii + ph 



(5.6) 



(5.7) 



and it allows us to better understand the difference between the new formulation (| 1 . 5|) and the classical 
one (jl.ip proposed by Coleman [5] and Truesdell and Noll [5^1151] . 

If, on the contrary, we want to keep the classical point of view, we must modify (jl.ll) by introducing 
a suitable extra flux N and assume (|1.7p . as suggested by Dunn and Serrin [5]. In fact, taking into 
account (|5.3I) - (|5.6|) and (|1.3I) . equation (|5.7I) assumes the form 



d (\ 



P dt e = -dt{2 PU 



pfii - V • N' + ph = T • Vii - V ■ (Tu + N') - V ■ q + pr, 



which coincides with (|1.7[) by choosing 

N = (Tu + N') = -a A u Vu. 

Let us now consider a thermoelastic plate with memory characterized by the following constitutive 
equation 

p f oo 

C'(s) Au t (s)ds + C Au(t) +cV9(t), (5.8) 



T(i) = -V 



where the scalar kernel C is a suitable smooth function and m*(s) = u(t — s). 

Reasoning in the same way as in the previous model, we evaluate the terms in the right-hand side 
of (|5.3p making use of (|5.8p and obtain 



d fl 



dt\2 



-pu 2 {t) 



C'(s) Au t (s)ds - c9(t) 



Au(t) + ^j t [Au(t)f 



j-V J^C'{s)Au t {s)ds + C Au{t) 
C'(s) A u t (s)ds + C A u(t) - c9(t) 



-cV6»(<)| u(t) 
Vu(t)) +pf(t)u{t). 



Hence, it follows that the internal and external mechanical powers are given by 

Co d 



KM) 

where 



C"(a) Au t (s)ds -c6(t) 



A «(t) + ^fj t [Au(t)Y , VIM = -V • N'(t) + p/(t)u(t) 



N'(t) 



{ V [P 



(s) Au 4 (s)ds + C Au(i) 



c\70{t) \ u(t) 



| Ml 



C'(s) A u t (s)ds + C A u(t) - c9(t) 



IJo 



Vii(t). 



Now, using the formulation (jl.5p of the First Law, we could define the internal heat power as follows 



Ph(t) - P -e(t) - 



C'(s) Au\s)ds -c0{t) 



a " w "tI [AuW]2 ' 



Finally, by choosing 



N(t) = - 



C"(s) A u t (s)ds + C A u(t) - cd(t) 



we can use also (JTTTJ) as local form of the First Law. 
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Other interesting examples, for which the classical representation of the power, used in the First 
Law, does not coincide with the one following from a correct balance between the internal and external 
powers, meet in electromagnetism, in the study of phase transitions and also in other physical systems, 
as those describing processes of phase separation, for example, by means of Cahn-Hilliard's equation 
0- 

Firstly, we consider a rigid electromagnetic system for which the classical expression for the First 
Law is given by 

— e = BH + DE + JE-Vq + r, (5.9) 
dt 

where E, H, D, B, J are the electric and magnetic fields, the electric displacement, the magnetic 
induction and the current density, respectively. 

Equation (15.91) together with (jl.3l) leads to the classical definition of internal electromagnetic power 

7^=BH + DE + JE (5.10) 

and thermal power 

h=— e-BH-DE-JE. (5.11) 
dt 

These relations hold for a large class of electromagnetic materials, such as classical dielectrics, conduc- 
tors characterized by the Ohm Law, materials with memory, etc. In analogy with the mechanical 
systems, these materials are still called simple. However there exist meaningful examples of fairly com- 
mon materials, for which (|5.9p must be modified by changing the internal electromagnetic power V\i 
or by introducing an extra flux as follows 

— e = B-H + D-E + J-E-V-N-V-q + r. (5.12) 
dt 

To this purpose we recall that, as already shown in |12j . in presence of superconducting materials it is 
necessary to use (|5.12l) with an extra flux proportional to the superconducting current. 

Here we consider a dielectric with quadrupoles described by the following constitutive equations 

D = e E - £1 A E - £ 2 V [V ■ E] , 

B = /iH, (5.13) 
J = 0, 

where £0, £1, £2 are suitable constants characteristic of the dielectric. 

The evolution of the electromagnetic field is governed by the well known Maxwell equations which, 
in absence of currents, are 

D = V x H, 
B = -V x E. 

An immediate consequence of the previous equations is the balance law for the electromagnetic power, 
known as Poynting theorem, 

D-E + B- H = VxH-E-VxE-H = -V-(ExH) 

which, taking into account the constitutive equations (|5 . 13[) . becomes 

E x H £1 VE E - £ 2 V -EE. 

Hence the expressions of the internal and external electromagnetic powers are 

Vh = ~ [ M |H| 2 + £ |E| 2 + £r|VE| 2 + £ 2 |V • E| 2 ] , 
E x H - £iVEE - £ 2 V • EE . 



-j [^IHp+folEp+^IVEp+^lV-El 2 ] =-V- 
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Consequently, the local form p.l[) of the First Law, with V l el instead of V^, becomes 



d Id 
dl e ~ 2~dl 



A*H 2 + e E 2 + e x (VE) 2 + e 2 (V • E) 1 



+ h 



and it allows us to give a properly definition of the thermal power h which does not agree with the 
classical definition (|5.11[) . 

Finally, it is easy to show that the relation (|5.12l) is satisfied by choosing 

N = -eiVEE - e 2 V • EE. 

Therefore, this example confirms that the local form (|3.1[) of the First Law (with V l el instead of V; n ) 
has a universal character; on the contrary, (|5.11[) must be adapted in function of the material under 
consideration. 

Let us now consider an example related to the phase transition in a binary mixture characterized 
by the Cahn-Hilliard equation 

c= V- [M(c)V/x(c)], ( 5 - 14 ) 

where c G [—1,1] denotes the relative concentration of one of the two components, the mobility M 
represents the density with which the two phases are mixed and fi is the chemical potential. 
The mobility M is a nonnegative function and the chemical potential [i is given by 

(i{c) = -jAc+6 F'{c)+6G'{c), (5.15) 

where 7 is a positive constant, #0 is the temperature of the transition, while G and F are the usual 
polynomials assumed for a second order transition, i. e. 

G(c)=/3y, F(c)=/3(j-£y (3>0. (5.16) 

Therefore, from (|5.14[) ~ (15.16l) . we obtain 

c = V ■ [M(c)V(-7 A c + e F'(c) + 6G'(c))] . (5.17) 
By multiplying (|5.17D by the potential /i, we have the equation for the power 

c>(c) = V • [M{c)^{c)V^{c)\ - M(c) [V^(c)] 2 , 
which, by virtue of (|5.16p . can be rewritten as 

e F(c) + 6G{c) + |— (Vc) 2 + M(c) [V^(c)] 2 = V • [ 7 cVc + M (c) //(c) V fJ,(c)\ . 
This relation allows us to define the internal and external chemical powers of the mixture 

VI = 9 F(c) + 06(c) + || (Vc) 2 + M(c) [V^cf , (5 18) 

V e c = V • pycVc + M(c)/z(c)V>(c)] • 

If we put (|5.18[) into the heat equation 

le-T*— V.q + r. 

we get 

= tfi(c) +7V • (cVc) + M(c) [Vfi(c)} 2 - V • q + r. (5.19) 

It follows that if we want to use the usual expression c/i(c) for the chemical power, even in this case, 
it is necessary to introduce the extra flux 

N = -7cVc. 

3 Among the wide literature on the Cahn-Hilliard equation we recall, for example, |28] and |19| that set 
precedents. 
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We conclude this subsection by observing that if we take for the internal energy the constitutive 
equation 

e(8, c, Vc) = 6 F(c) + 1 (Vc) 2 + e(0), 
then equation (|5.19[) assumes the following more usual form 

^(0)6 = 6G(c) + M(c) [V^(c)] 2 - V • q + r. 

A similar argument can be also applied to the binary mixture of two incompressible fluids studied 
by Lowengrub and Truskinovsky in |24) . 

5.1.2 Second Law of Thermodynamics 

In this subsection we consider only non-simple materials in the heat flux, in other words, materials for 
which the constitutive equation for the heat flux depends also on the gradients of order higher than 
the first in the temperature. 

For second grade materials, taking in account Definition 15. 11 the internal virtual entropy action 
becomes 



4„ (a, P, 7 t ) = ph(a, P) i - qi (a, P) ■ V 



1 



q 2 (a, P)-V< a > 



1 



Inequalities (jl.2p and (|1.10[) or, equivalently, 
d h 1 „ n 



for simple materials, 



p— r\ > p — | — -q ■ V0 — V • 3»' for non-simple materials, 
dt 8 8 2 



(5.20) 

(5.21) 
(5.22) 



show that the Second Law assumes different expressions for simple or non-simple materials. 

Reasoning in same way as for the First Law, we suggest to claim the Second Principle in a unique 
general form trough the inequality (j!.17[) . 

For simple materials the internal entropy action A\ n is always given by the differential form of the 
right-hand side of (|5.21j) . while the right-hand side of (|5.22D cannot be taken as expression of A l en for 
non-simple materials, because of the presence of the external flux V • 

We now present an example [T] that helps to better understand the motivations of our suggestion. 
Let us consider a rigid heat conductor, whose constitutive equation is given by 

q(t) = qi(0(t),§*) - V- q 2 (0(t),g\ Vtf) (5.23) 

where qi is a vector, q 2 is a second order tensor and g* is the integrated history of V0, that is 



r(X)d\. 



The functional (|5.23l) is a typical constitutive equation of a second grade material; therefore, taking 
into account (|5.20p . the internal entropy action becomes 



A' 



1 



Oh + qi • V0 + q 2 • V (2) 6» - -q 2 V6> • V0 



Consequently, the First and Second Thermodynamic Laws have the following expressions! 
d 



dt 



= h. 



d 1 
dt^ ¥ 



6h + qi • V0 + q 2 • V (2) - -q 2 V0 • V0 



(5.24) 
(5.25) 



4 We can assume p = 1 without loss of generality. 
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We observe that, choosing 



the inequality (|5.22|) coincides with (|5.25[) . 

Finally, we assume the following constitutive equations for the internal energy and heat flux 



e(0(t),t,V&) =e(0(<))+6(g t ,Vg*), 

r+oo 

qi(0(i),g') =-0(t)/ 1^ («)!'(«)*», 

./o , 

r+oo 

q 2 (e(t),g t ,Vg t ) = -e(*) / K' 2 (s)Vg t (s)ds, 

Jo 



with e positive and homogeneous function of 9, K\ and K% smooth positive functions. 
Then the internal entropy action becomes 



4„W 



i 



9{t) 



r+oo r-t-oo 

h(t) - / K[(s)g t {s)ds ■ V0(i) - / K' 2 (s)Vg t (s)ds ■ V {2) 9(t) 
Jo Jo 

K! 2 (s)Vg t (s)dsV9(t) ■ V0(t) 



^(*) Jo 

and the free energy potential (|4.6[) can be expressed as a sum of two terms 

^(0(t),g*, Vg*) - iM0(t)) + Ms\ vg 4 ), 

where 



and the functional "02 satisfies the inequality 

d r+ °° 



ld_ 

t dr 



e(r)dT 



dt 



p+oc r+oo 

^(gSVg*)) < / K' 1 (s)g t (s)ds-V9(t) + K' 2 (s)Vg t (s)ds-V^9(t) 
Jo Jo 



5.2 Other non-simple materials 

It is worth to observe that there exist material systems, whose constitutive equations present higher 
order gradients, for which the internal mechanical virtual power and internal virtual entropy action 
are not covered by the expressions (|4.ip and (|4.2|) . 

For this purpose, let us consider the model proposed by Guyer and Krumhansl [2T] for the heat 
flux and recently studied in [4 ,32]. Since the constitutive law for heat flux is given by 



q + — q = -c v V9 + t n [Aq + 2V (V • q)] , 



(5.26) 



it is easy to show that for this material the internal virtual entropy action cannot be represented in 
the form (|4.2j) . Anyway, we can still define the internal entropy action. In fact, if we choose 



c o ^ n 

Cv = ~g~2> c ° > °> 



and multiply (I5.26P by q, we obtain 



J_ q . V = -l 

9 2 c 



1 d 
2~dt 



— q 2 + tat (Vq) 2 + 2 (V • q) 2 
tr 



Co 



V- (Vqq + 2Vqq). (5.27) 
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The heat equation for a rigid conductor together with (|5 . 2T[) yields 



Id 1 d 2 
8di e ~ 2co"dt q 



1 2 T N 

q' 



(Vq) 2 + 2 (V • q)' 



-V- (2) +^-^V-(Vqq + 2V-qq). 



Hence, the expressions of the internal and external entropy actions are given by 



A' 
A' 



Id 1 d 2 
di 6 ~ 2co"dt q 

r 



1 



q 



(v q ) 2 + 2 (v • q y 



where 



*d = — (Vqq + 2V-qq) 

CO 



and the Second Law of Thermodynamic (I1.17[) becomes 



1 d 



1 d 



-nV>^-re-——a 



dt 



9 dt 2c dt 



CqTr Cq 



(Vq) 2 + 2 (V • q) 



(5.28) 



The previous inequality imposes restrictions on the constitutive equation (|5.26p . In fact, if we choose 
the following expression for the entropy 

" 1 de , 1 2 
T OT 2cq 



then from (|5.28l) we obtain that the inequality 



1 



TN 



(Vq) 2 + 2(V-q) 



> 



coTr c 

must be satisfied for any held q and therefore the constitutive coefficients tr and tn must be positive. 
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